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We report magnetoresistance measurements of ramp type superconductor-normal-superconductor 
(SNS) junctions. The junctions consist of underdoped Y BaiCuzOy (YBCO) electrodes separated 
by a barrier of YBa2Cu2.6Gao.4Oy. We observe a large positive magnetoresistance, linear in the 
field. We suggest that this unusual magnetoresistance originates in the field dependence of the 
proximity effect. Our results indicate that in underdoped YBCO/N/YBCO SNS structures, the 
proximity effect does not exhibit the anomalously long range found in optimally doped YBCO 
structures. From our data we obtain the diffusion coefficient and relaxation time of quasiparticles 
in underdoped YBCO. 
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PACS numbers: 74.45.+C , 74.25.Ha 

INTRODUCTION 

In the usual description of the proximity effect, when a 
superconductor (S) is brought into contact with a normal 
conductor (N), the order parameter (OP) in the super- 
conductor is depressed near the interface and supercon- 
ductivity is induced in N. The pair amplitude induced in 
N decays on a length scale K~ x from the interface, called 
the decay lengthy], 0,y]. In SNS junctions in which S is 
an optimally doped HTSC and N belongs to the same ma- 
terial family, but is doped to be non superconducting, the 
decay of the pair amplitude in N typically takes place over 
a rather long distance of tens of nm,0,|{|. In contrast, if 
both S and N are underdoped cuprates, the pair ampli- 
tude in N seems to decay over a much shorter distance, 
on the order of a few nm. We have observed this effect in 
underdoped YBCO based junctions0,Q having a barrier 
made of YBa2Cu2.55Feo.5Oy or a YBa2Cu2.6Gao.4Oy. 
In SNS junctions having a barrier much thicker than 
the decay length, Cooper pairs cannot tunnel through 
and the junctions exhibit a finite resistance at all tem- 
peratures. Roughly speaking, superconductivity in N is 
induced near the two SN interfaces, while a section of 
length t in the middle of the barrier remains normal. 
This is the type of junction studied in the present work. 

We are not aware of previous investigations of the prox- 
imity effect in HTSC under a magnetic field. When a 
magnetic field is applied, superconductivity is reduced 
and penetrates less into the normal conductor. As a 
result, the proximity effect is field dependent 0. If 
the superconductivity in the barrier is weakened, the 
length of the normal section in the junction should in- 
crease, and with it the junction's finite resistance. As 
a result, a positive magnetoresistance (MR, defined as 
MR = R(H) - R(0)) should be observed. We indeed 
observed such MR, linear in the field. An attempt to ex- 
plain this unusual field dependence is the subject of this 
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paper. 

Besides the field dependence of the proximity effect, 
there are several additional mechanisms which could 
contribute to the MR. These include flux flow in the 
superconducting electrodes 0, EJ, normal MR of the 
barrier material which is caused by bending of electron 
trajectories [TH. field dependent hopping in the barrier 
[l2l Il3| , and resonant tunneling between the electrodes 
across the barrier [l^. In the following we show that 
the contribution of all these processes to the observed 
MR is insignificant and we attribute it primarily to a 
field dependent proximity effect in the barrier. 

EXPERIMENTAL 

The junctions used in the present study are thin film 
based ramp junctions of the type that was previously 
used in our work The junctions consist of two un- 
derdoped superconducting YBCO electrodes separated 
by a thin layer of Ga-doped YBCO barrier. Ga has 
no magnetic properties. The transport current flows 
in the a-b plane through the Ga-doped YBCO layer. 
The multi-step process of junction preparation by laser 
ablation was described previously 6] . Briefly, we first 
deposit a 100 nm thick c-axis oriented epitaxial YBCO 
layer onto a (100) SrTi0 3 (STO) substrate. This base 
electrode is then capped by a thick insulating layer of 
STO. Patterning is done by Ar ion milling to create 
shallow angle ramps along a main crystallographic 
direction in the a-b plane. In a second deposition step, 
the barrier layer, the YBCO cover electrode and the Au 
electrical contacts are deposited, and then patterned to 
form the final junctions layout. This produces several 
junctions with 5/im width on the wafer. Four terminal 
resistance measurements of the junctions were carried 
out as a function of temperature and magnetic field 
of up to 8 Tesla. The field was perpendicular to the 
transport current, which in our geometry flows in the 
a-b plane of the films. 
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FIG. 1: (Color online) Resistance versus temperature of six 
junctions with 21nm thick Ga-doped YBCO barrier. In the 
normal state, the different resistances of the junctions are 
due to different lengths of the YBCO leads. The inset shows 
the low temperature resistance of the junctions where both 
electrodes are superconducting. 
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FIG. 2: (Color online) Resistivity versus temperature of 
lOOnm thick film of the Ga-doped YBCO material. Square 
symbols are measured at zero field and the circles are mea- 
sured with 6 Tesla field applied perpendicular to the film. 
Note that the MR of the film is negative, in contrast to the 
positive MR of our junctions. 



RESULTS AND DISCUSSION 

Resistance versus temperature (RT) measurements 
of six junctions on the wafer are shown in Fig. 1. 
In the normal region, the difference in the resistance 
of the junctions is due to the different lengths of the 
YBCO leads. One observes two distinct superconducting 
transitions with T c onset of 35K and 53K which are 



attributed to each one of the electrodes. In the oxygen 
annealing process of underdoped YBCO, the oxygen 
concentration is kept low and the duration of the anneal- 
ing is relatively short. Consequently, the base electrode 
which is covered by a thick layer of STO, absorbs less 
oxygen and its transition temperature is lower. Below 
about 30K, both electrodes are superconducting and the 
inset of Fig.l shows the low temperature resistance of 
the junctions, which is due to the barrier. Qualitatively 
similar behaviour was observed in edge junctions made of 
underdoped YBCO separated by a YBa2Cu2.55Feo.5Oy 
barrier ja, 0- The scatter of the values between different 
junctions is typical of our junction preparation process, 
and is probably due to nonuniformities in the local Ga 
concentration, and to variations in the transparency of 
the interfaces, most probably resulting from damage 
created by the ion milling of the ramp. The transparency 
of our junctions can be estimated from measurements of 
the critical current described below, which indicate that 
the transparency is low. The temperature dependence 
of the junction's resistance is typically weaker than that 
of the parent material in the form of a film, shown in 
Fig. 2. At low temperatures, the absolute resistivity 
of the junctions is also much smaller than that of a 
YBa2Cu2.6GaQ.4Oy film. One possible interpretation is 
that the thickness of the barrier (21 nm in this work) is 
in the range where the material is mesoscopic. Under 
these conditions, the temperature dependence of the 
resistance is expected to be much weaker than that of 
a macroscopic filrn^^]. The differences of the absolute 
resistivities between different junctions may perhaps 
result from different interface transparency, which also 
affect the conductance of the device in the mesoscopic 
regime 15]. 

Our main experimental result is shown in Fig. 3, 
where the measured magnetoresistance MR at low T is 
plotted as a function of magnetic field normal to the 
wafer. All junctions showed similar behavior. Detailed 
measurements were done on three of the six junctions 
on the wafer. One can see that all three junctions show 
a large positive MR which is linear in the applied field. 
The MR typically reaches a value of ~ 20S1 at 8 Tesla, 
which is larger than the resistance at H=0 by tens of 
percent. 

We consider possible sources for the MR in our 
junctions. MR originating in the two YBCO electrodes 
below the superconducting transition temperature can 
result for instance from flux motion. This contribution 
would be linear in the field. In order to estimate the 
size of this contribution, we performed low temperature 
MR measurements on bare YBCO microbridges. At 
temperatures close to T c , flux flow was indeed observed 
(see Appendix ). However, at low temperatures where 
the junctions of Fig. 3 were measured, no measurable 
MR was observed in the thin film YBCO microbridges. 
Therefore, flux flow in the YBCO electrodes does not 
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at most -851 (5% of 160J1, as seen in the inset of Fig. 
1). However, since the sign of the MR of the barrier 
material itself is negative, the net (positive) MR of the 
junctions should be even larger than shown in Fig. 3. 
Consequently, the properties of the barrier material on 
its own cannot explain the observed MR of the junctions. 
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The above mentioned control experiments clearly show 
that the MR of our junctions does not originate from 
the superconducting electrodes nor from the normal 
properties of the barrier material. The net MR which we 
see has a magnitude characteristic of the transition of 
part of the barrier from a superconducting to a normal 
state. We therefore examine whether the MR could 
originate from the depression of superconductivity near 
the SN interface of the junction. 



FIG. 3: (Color online) Magnetoresistance versus field of three 
of the junctions of Fig. 1 at 2K, and of a "short" junction 
at 4K. The "short" resistance is of about 0.4Q at 8 Tesla 
which is almost two orders of magnitude smaller than the 
corresponding MR of the other junctions with a barrier. The 
dashed line is a guide to the eye. 



contribute to the MR. We also measured the MR 
of junctions prepared in the same way, but without 
the barrier layer. These junctions are referred to as 
"shorts". As shown in Fig. 3, under similar bias 
currents and fields, the "shorts" did exhibit a small 
MR of about 0.417 at 8 Tesla. The "shorts" show a 
finite MR since the interface between the two YBCO 
electrodes is always imperfect. The MR of the "shorts" 
is smaller than the MR in the junctions by almost two 
orders of magnitude. The interface resistance cannot 
be directly measured. What can be measured is the 
critical current density. Typically, the critical current 
density at low temperature of a 60K YBCO "short" 
is one order of magnitude smaller than that of a film. 
This implies that the transparency of our junctions is 
low. To summarize this section, the above series of 
control experiments show that MR in the electrodes is 
not the source of the large MR observed in our junctions. 

A second potential source for the observed MR in 
Fig. 3 could be the barrier material itself. We therefore 
measured the MR of the Ga-doped YBCO. Specifically, 
we measured the resistance versus temperature of 
microbridges patterned in a thin film of this material 
annealled under the same conditions as the junctions 
in Fig. 3. Fig. 2 shows the resistivity of these bridges 
with and without magnetic field. The barrier material 
exhibits a clear negative MR of ~5% at 2K. The sign 
of this MR is opposite to that of the junctions which 
show a large positive MR. At low temperatures, where 
the MR of the Ga-doped films is largest, the MR 
contributed by the barrier in the junctions would be 



Before going into a more detailed analysis, we note 
that our barrier is a mesoscopic section of a Mott 
insulator (MI), with the conductance of the material 
in bulk form showing variable range hopping |16| . Its 
low temperature resistivity, 0.8J7cm, is about 3 orders 
of magnitude larger than the maximum resistivity of 
metals (Mott-Ioffe-Regel limit [2^). Strictly speaking, 
our junctions are S/MI/S junctions. So, the application 
of the usual theoretical description of the proximity 
effect to our junctions is not a priori justified, since 
both the de-Gennes and Usadel equations are valid 
only for dirty metals. However, it is an experimental 
fact that when an MI with resistivity p < 1 flcm is in 
good electrical contact with a superconductor it behaves 
similarly to a metal 123, l24l l2q . The question which 
particular model to use is therefore a matter of choice. 
In the limit of small induced pair amplitude in N, which 
applies to our low transparency junctions, the de-Gcnncs 
and Usadel approach give the same result. Since the 
de-Gennes approach was traditionally employed in 
all previous and current work on HTSC proximity 
structures^ |U|2(|, we prefer to follow this route. In any 
case, the analysis presented below is nevertheless useful, 
in terms of assigning values to physical quantities such 
as the decay length which can then be inter-compared 
between different experiments. 

We first discuss the MR on the S side of the SN 
interface. In this region, the order parameter is reduced, 
superconductivity is depressed, pinning is weakened 
and flux flow could occur despite the low temperature. 
We now estimate the upper limit on the contribution 
of this effect to the MR. The low temperature normal 
state resistivity of YBCO, extrapolated from the linear 
part of the RT plot above the transition, is about 
10 -4 fic77i. An upper limit on the volume near the 
interface in which superconductivity is weakened is 
10£ xi~ 200Ax0.5/im 2 , where A is the junction cross 
section [2lJ • The normal state resistance of this region is 
very low, less than 0.112. Since the flux flow resistance is 
a fraction of the normal state resistance, it follows that 
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the MR in the S side close to the interface is negligible. 

Turning now to the N side of the interface, the 
resistivity of the barrier material is quite high, 0.8f2cm 
at 2K. A rough estimate done assuming Ohm's law 
in the barrier indicates that a 1 nm thick slice of the 
barrier has a resistance of R ~ 160. This value is similar 
to the total MR seen in Fig. 3. In the following, we 
propose that the observed MR is caused by changes 
in the effective penetration of superconductivity into 
the barrier. In other words, when a magnetic field is 
applied, the magnitude of the pair amplitude induced in 
the barrier is decreased and £, the effective length of the 
barrier which remains normal, increases thus increasing 
the resistance of the junction. 

The magnetic fields used in the present study are small 
compared with H C 2 of the 60K YBCO phase which is 50T 
|2fij . Thus changes in the minigap A due to the applied 
field are also small but not negligible. The value of A on 
the S side near the interface is proportional to T c which 
itself depends on the magnetic field due to pair breaking 
according to 0, 0] . 



Proximity 




FIG. 4: A schematic diagram of the junction and the spatial 
profile of A(x). 1(H) is the length of the resistive region of 
the junction. The shaded area shows the region in which 
superconductivity is weakened on both sides of the interface 
due to the proximity effect. 
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where T c is the critical temperature under applied field 
and T c q is the critical temperature at zero field. ^ is the 
di-gamma function defined as ^>(x) = T'(x)/T(x), and 
a is the pair breaking parameter. For a thin film under 
perpendicular applied field a = DseH/c, where D$ is 
the diffusion coefficient in the superconductor. Since the 
highest magnetic field we used is small compared with 
H C 2, pair breaking is small and (a/27rfcsT c ) is a small 
parameter. In this limit, Eq. (1) reduces to 
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where A s and A l N are the values of the minigap at the 
S and N sides of the SN interface. N$ and Njy are the 
normal state density of states (DOS) on the S and N 
sides of the interface, respectively. Finally, Vs and Vjv 
are the electron-electron interaction on the S and N sides. 
Assuming the DOS and the electron-electron interaction 
are field independent we obtain 



ks{T c0 - T c ) 



(2) 



From our RT measurements under different fields we find 
the values of T c0 and T c (65K at OT and 55K at 7T, 
respectively). We can thus calculate the value of a which 
is ~1 meV at 7 Tesla. Therefore, a/27rfcsT c ~ 1/35, and 
this justifies the use of Eq. (2). By assuming a linear 
scaling between A and T c (2A = /3ksT c , with (3 being 
a constant of about 5), we estimate that under a field 
of 7T the magnitude of A decreases by about 15%. The 
suppression of A can be therefore written as 
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where e = NsVs/NnVn is a field independent constant 
and we define Sn = Aat(O) — An(H). 5 l N which is 
the value at the interface, is also a small parameter as 
VAivW « 1. 



Turning to the N side now, under a magnetic field H 
applied in the c-direction, the spatial dependence of A is 
given by the linearized Ginzburg-Landau (GL) equation 



5 S = A s (0) - A S (H) 



iraf3 7T/3 D$eH 
~ ~ ~8 c 



(3) 



where 5s is small compared with Ag(0). The spatial de- 
pendence of the A in an SNS junction is shown schemat- 
ically in Fig. 4. The value of the A on both side of 
the interface are related through the standard boundary 
condition 
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(x - x) 2 A N + K 2 A N = (6) 



where xq and K are constants. In our experiment, xq — x 
is limited by 10 nm which is half the thickness of our 
junction, the field H is less than 8T, and K is on the order 
of a few nm. Using these parameters, we estimate that 
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the upper limit of the second term in Eq. (8) is about 
two orders of magnitude smaller than the last term. In 
this limit, the solution of Eq. (8) for A in N exhibits an 
exponential decay with distance An (x) — A N exp(—Kx). 
In the dirty limit [2| , K is given by 
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N 



2nk B T 
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In this limit, where the cyclotron radius in the magnetic 
field is much larger than the mean free path, Dn is field 
independent and thus K does not depend on field. How- 
ever, the value of A at the interface A l N , is field depen- 
dent because it is pinned to the value of the A on the S 
side at the interface A l s through Eq. (5). The pair am- 
plitude induced in the barrier is effectively depressed to 
zero by thermal fluctuations at some distance from the 
interface, and from that distance onwards the material 
has a finite resistance. The natural way to determine 
this distance is through the condition that the extrapo- 
lated magnitude of A there is of the order of k B T . This 
length, which we denote by A, depends on the field as 



(8) 



where X(H) is the effective penetration depth of super- 
conductivity into N when a magnetic field is applied. 
Dividing the Eq. (9) for X(H) by the equation for 
X(H = 0) we find 
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N is a small parameter X{H) 



X(0) 



—S N /KA N (0). Referring to the schematic model shown 
in Fig. 4, the field dependent resistance of the barrier is 
R = p£(H)/A, where 1{H) = L - 2X{H) is the length 
inside the barrier which is normal. Using Eq. (5) and 
the relation 2As = f3ksT c , the magnetoresistance comes 
out as 



MR = R(H)-R(0) = -2^(X(H)-X(0)) (11) 
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We therefore see that the MR is linear in H, in agreement 
with the observed behavior in Fig. 2. 

A rough estimate of the decay length (1/K) in 
the underdoped barrier at low temperature can be 



attempted using the resistivity of the barrier, 0.8f2cm 
and the typical resistance of the junctions ~ 100S1. 
Using these values, we estimate the length of the barrier 
which remains normal £(H — 0) as 6nm. Taking the 
thickness of the barrier of 21nm, and assuming that 
the pair amplitude decays to zero over 3 times the 
decay length (1/K), we obtain a value for 1/K of about 
2.5nm. 1/K can also be calculated using Eq. (7), where 
Dn — \(-nVfn- The mean free path in the barrier 
can be estimated as the distance between nearest Ga 
atoms £n ~ 5A and the Fermi velocity in the barrier 
vfn = 1-2 x 10 7 cm/ sec was measured in a previous 
study[3. This yields 1/K ~3.5nm. It appears that 
both methods of estimating 1/K give values which are 
consistent. We note that the decay length estimated in 
underdoped SNS structure comes out much smaller than 
in optimally doped ones 



Using our data we estimate the diffusion coefficient 
Ds and the relaxation time ts of underdoped YBCO. 
Taking an average T c of ^45K, and an average slope 
in Fig. 3 of MR/H = 2.7Q/Tesla, Eq. (11) yields 
Ds ~ 1cm 2 / sec. This is consistent with an independent 
estimate that can be extracted from Eq. (2) and from 
the relation between a and the diffusion coefficient Ds 
which yields ~ 1.7cm 2 / sec. The relaxation time, ts is 
extracted from the usual relation that connects it with 
the diffusion coefficient Ds — ^Vp S Ts where Vfs is the 
Fermi velocity of quasiparticles in the superconductor 
~ 2 x 10 7 cm/ sec 27]. Under these assumptions ts 
for YBCO is ~ 25/s. The value found for ts is of 
the same order of magnitude as the recent results 
of Gedik et al. who obtained ts ~ 100 fs |27j while 
our value oi Ds is smaller by than theirs, Ds ~ 20 cm 2 /s. 



For completeness, we mention that Abrikosov has 
predicted another mechanism for linear MR versus H 
in superconductors |l4j |. He assumed a field dependent 
resonant tunneling which yields MR linear in H at very 
high magnetic fields, when only a few Landau levels 
are filled. When the field is reduced and the number 
of filled Landau levels increases the field dependence of 
the MR changes into a quadratic one. This model could 
in principle, explain the observed linear behavior of our 
MR results. However, peaks in the density of states due 
to Landau levels are absent in the dynamic resistance 
spectra of our junctions. Moreover, the fields used in 
our experiment are not high enough to reach the regime 
where a low number of Landau levels arc filled. Hence, 
if this model was applicable to our junctions, we should 
have observed a quadratic dependence of the MR on 
field, which is not the case. 

CONCLUSIONS 

We investigated the resistance of SNS structures based 
on underdoped YBCO with non-magnetic Ga-doped 
YBCO barrier as a function magnetic field. We discov- 
ered a linear increase of the resistance with the field 
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FIG. 5: (Color online) H C 2 versus temperature of optimally 
doped YBCO film. H C 2 was extracted from the slope of 
the linear part of the MR (inset) using the Bardeen-Stephen 
model. Our data (solid squares) can be compared with pre- 
vious measurements by Kunchur et al. and Ossandon rs 



MR. An extensive series of control experiments indicates 
that this field dependence does not result from flux flow, 
which would be the obvious mechanism of MR in a 
superconductor. A simplified analysis indicates that the 
effect may well be explained by field dependent prox- 
imity effect in the barrier. This explanation produces a 
reasonable estimate of the diffusion coefficient and the 
relaxation time in YBCO. Furthermore, our estimates 
indicate that in underdoped YBCO SNS structures, the 
superconductivity induced inside the barrier through the 
proximity effect has a short (~2-3 nm) range, unlike the 



long range proximity effect observed in optimally doped 
YBCO structures. 
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A. APPENDIX 

For the sake of comparison with previous work, we 
also measured the MR of optimally doped YBCO films 
at temperatures close to T c as shown in the inset of 
Fig. 5. In this case, the resistance showed a region 
linear with applied field. In the Bardeen-Stephen model 
El 113 5 the resistance resulting from flux flow is given 
by R f inx flow = (H/H c2 {T)) x R N (T) where H is the 
applied magnetic field and Rn(T) is the normal state 
resistance at temperature T, extrapolated from the RT 
plot close to T c . Using this model, we extracted the 
temperature dependence of H C 2 near T c . Our results 
show good agreement with previous measurements by 
Kunchur et al. and Ossandon et al. [2^, [3(j > which are 
also plotted in Fig. 5. At temperatures much lower 
than T c however, no measurable MR in the YBCO film 
was observed. Therefore, at low temperatures where 
the junctions of Fig. 3 were measured, flux flow in the 
YBCO electrodes does not contribute to the MR. This 
conclusion holds independent of the oxygen doping level 
of the YBCO. 
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